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TECHNICAL NOTE 2j8j 


AIRFOIL PROFILES FOR MINIMUM PRESSURE DRAG AT SUPERSONIC 
VELOCITIES - APPLICATION OF SHOCK-EXPANSION THEORY, 
INCLUDING CONSIDERATION OF HYPERSONIC RANGE 
By Dean R. Chapman 


SUMMARY 


A theoretical investigation is made of airfoil profiles at supersonic 
velocities to determine the shapes having minimum pressure drag at zero 
lift for various given auxiliary conditions. Shock-expansion theory is 
employed, thereby extending the applicability of the results through the 
hypersonic range. Curves are presented for Mach numbers of 1.5, 2, 3, 

L, 6, 8, and » which enable the Bhape and the drag of an opt imum profile 
to be determined readily if the base pressure is known from experiments. 
Examples are presented of optimum profiles determined with the aid of 
experimental base pressure data. Variations in profile shape are inves- 
tigated to provide information on the degree to which deviations in 
shape from the optimum can be made without resulting in a significant 
drag increase. 

A comparison of opt imum profiles determined by the shock- expans ion 
method of this report with corresponding profiles determined by the 
linear ized-theory method of a previous report shows only small differ- 
ences in shape at Mach numbers up to infinity even though the linearized 
theory at high supersonic Mach numbers breaks down completely insofar as 
the drag of the profile is concerned. The experimentally observed 
dependence of base pressure on trailing-edge thickness is found to have 
a significant effect on the shape and drag of optimum profiles of sma 3 ~l 
thickness ratio. Curves are presented which show that for thin airfoils 
the use of a trailing-edge thickness considerably greater than the theo- 
retical optimum can result in an excessive drag penalty at moderate 
supersonic Mach numbers, though not at hypersonic Mach numbers. 


INTRODUCTION 


In 1933 Saenger observed that for the extreme case of flow at 
infinite Mach number an airfoil designed to have minimum pressure drag 



would have its maximum thickness at the trailing edge. (See refers 
ence 1.) A related result can be inferred from, the numerical calcula- 
tions of Ivey (reference 2} which indicate that the drag of a 10-percent- 
thick-wedge airfoil at a Mach number of 8 is less than that of a double- 
wedge airfoil having the same thickness ratio. In both of these cases 
the desirability of employing a thick trailing edge in conjunction with 
a small surface slope may be attributed to the fact that at hypersonic 
Mach numbers the suction forces (forces due to pressures below ambient) 
are small compared to the positive pressure forces, even when the suction 
force corresponds to a vacuum. Recently, Smelt (reference 3 ) has dis- 
cussed this latter characteristic of hypersonic flow and its possible 
application to the determination of efficient airfoil shapes for use at 
very high Mach numbers. The investigations _gf Saenger, Ivey, and Smelt, 
however, do not provide general quantiiative information on the airfoil 
profile having minimum pressure drag in hypersonic flow because of two 
limitations of their analyses: Airfoils having a trailing-edge thickness 

less than the maximum airfoil thickness were. not considered, and the 
airfoil structural characteristics were notrconsidered (comparisons were 
made on the basis of a given airfoil thickness). 

At low and moderate supersonic Mach numbers the suction forces on 
an airfoil can amount to several times the positive pressure forces, 
particularly if an excessively thick trailing edge is employed. As a 
result, the optimum trailing-edge thickness t in this lower Mach number 
range depends to a great extent on the base pressure. By presuming that 
the base pressure is known from experiments, and that the airfoil profile 
must satisfy a given structural requirement (such as a given section 
modulus or a given section moment of inertia), a method of calculating 
the profile of minimum pressure drag at zero lift has been developed by 
the present writer in reference 4. Although the basic equations devel- 
oped in reference 4 for calculating such profiles are. applicable to 
higher-order theories, a detailed solution was given only for the case 
of linearized supersonic flow. 

Because of the well-known shortcomings of linearized theory, it was 
thought worthwhile to conduct an investigation parallel to that of 
reference 4 in which the shape and drag of optimum airfoils are determined 
from shock-expansion theory instead of linearized theory. Shock-expansion 
theory appears adequate for this purpose, particularly in view of the 
recent investigation of Eggers and Syvertson (reference 5 ) which indicates 
that shock-expansion calculations accurately determine surface pressures 
on thin airfoils in inviscid flow at Mach numbers from just above that 
for bow-wave attachment to infinity. From the viewpoint of the engineer 
who always has to make design compromises, it was thought desirable in 
the present study also to determine how much the optimum -profile shape 
can be altered, especially near the trailing edge, and still not increase 
the drag excessively. The purposes of the present investigation, there- 
fore, were (l) to develop a usable method for determining the shape and 
drag of optimum profiles in the Mach number range beyond that covered 
adequately by linearized theory (step-by-step details involved in applying 
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the method developed are' given, in an appendix), and (2) to determine 
curves shoving the rate at which the total pressure drag increases as 
the profile shape deviates from the optimum. 


NOMENCLATURE 


■t 


c 

c d 

C 

h 

H 

I 

T 

Jn 

l 

L 

M 

n 

P 

Pt 

P 

Re 

s 

t 

V 

x 


airfoil chord 
pressure drag coefficient 
constant depending on 7 
trailing-edge thickness 




dimensionless trailing-edge thickness ^ J 

given value of auxiliary integral 

I 


dimensionless value of I 


{■f [(t/2) a J 


jn-cr 




l(t/2)= 

normalizing factor for ?(Y), defined by equation (lo) 

length of surface of constant thickness 

dimensionless length of surface of constant thickness (— 

V s . 

Mach number 

parameter appearing in definition of I 
pressure 
total pressure 


pressure coefficient 


'p^ jPg\ 


Reynolds number 

chordwise distance from leading edge to first downstream position 
of max im um thickness 

maximum thickness of airfoil 

velocity 

chordwise distance from leading edge to point on airfoil surface 
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y 

X 

Y 

r 

x 

p 

8 

l(Y) 

t](Y) 


0 

1 

00 

b 

ca 


1 


ordinate of upper surface of airfoil 


dimensionless distance I — 

\ s 

dimensionless distance 

ratio of specific beats (1.40 for air) 

arbitrary constant 

mass density 

local angle of inclination of airfoil surface with respect to 
chord line (tan“ 1 y’ ) 

characteristic function determining optimum-profile shape, defined 
by equation ( 9 ) 

characteristic function ^defined by equation (12) 


Subscripts 


airfoil surface at leading edge 
airfoil surface at trailing edge 
free stream 

base, or trailing edge, 01 airfoil . 

circular-arc biconvex airfoil having sharp trailing edge 

Superscripts 


differentiation with respect to x 


ANALYSIS 


Solution for Arbitrary Structural Requirement 

As in reference 4> it is assumed throughout this analysis that the 
optimum airfoil has a sharp leading edge, a fixed. chord length, and is. 
set. at the zero-lift angle. It also is assumed that the flow is a 
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purely supersonic two-dimensional flow of an inviscid, nonconducting, 
perfect gas. Since the surface pressures on the top and bottom of an 
airfoil can be calculated independently in a supersonic flow, it follows 
that at zero lift the optimum profile will be symmetrical about the 
chord line. The mathematical problem formulated is to find the airfoil 
ordinate function y(x) which minimizes the pressure drag, 

2 f e h ' , , 

Cd = c Jo p y' ** - ^ C (!) 


for a given value of the auxiliary integral 


I 


_ 1 f c __£ L. 

C J o 


dx 


(2) 


By selecting various values of the parameters n and a, a wide variety 
of structural requirements can be represented for both thin-skin and 
solid-section structures. Same of the different structural criteria 
represented by equation (2) are: 


n 

a 

Structural criteria 

1 

0 

given torsional stiffness, or torsional strength, of 
thin-skin structure (given cross-section area) 

2 

0 

given bending stiffness of thin-skin structure 

3 

0 

given bending stiffness or given torsional stiffness of 
solid-section structure 

2 

1 

given bending strength of thin-skin structure 

3 

1 

given bending strength of solid-section structure 


Basic equations.- The equations which the optimum-airfoil ordinate 
function y(x) must satisfy can be obtained by considering an infinitesi- 
mal variation in profile shape 5y(x) that is arbitrary except for the 
requirement that 51 = 0. By also requiring that 5cd — 0 the following 
three equations result (see reference 4): 


s( p + y, IO ' +wn ' 1 ' 0 

(3) 

/dp \ 

+ Wvi 

W 

I = 2 T (T = i/(t/2) n -°} 

(5) 


In the derivation of these equations it was assumed that P .= p(y* , yo T , M^), 
but no particular functional form was assumed. Equation ( 3 ) is the differ- 
ential equation which the airfoil ordinate function y(x) must satisfy along 
the curved surfaces, for example, 0A and BC in figure 1. A first integral 
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of this differential equation, satisfying the condition y* = 0 at 
y = t/ 2 , can he obtained by multiplying both. sides by y* and integrating 

= \[(t/ 2 ) n - y 1 ] ( 3 a) 

Equations (3) and (3a) involving the arbitrary constant X do not apply 
to the straight midsection (AB in fig. l) along which the airfoil thick- 
ness is constant. Equation (1*.), termed an end condition, represents the 
relation which must be satisfied between the base pressure coeffi- 
cient Pb, the surface pressure coefficient P 1 Just upstream of the 
trailing edge, and the corresponding surface slope y 1 *. As will be 
seen this equation determines the optimum trailing-edge thickness. 
Equation (5) relates the optimum length of straight midsection l to 
the dimensionless value of the structural integral I. This latter 
equation ■ shows that the length of straight midsection is always zero 
when the auxiliary condition represents a given stiffness (cr = 0 ), but 
for the values of n considered is a sizable fraction of the chord 
when the auxiliary condition represents a given strength (cr = l). 


Solution for Y, H, and X .- From equation (3a) the constant X is 
readily evaluated in terms of y 0 *. There results 

k(t/ 2 ) n = yo , 2 (^p^ (3b) 

By employing the dimensionless variables Y = y/(t/2) and H = h/t, 
equations ( 3 a) and ( 3 b) yield 


y’^dP/dy* ) 

y 0 ,2 ( Sp /^y’)o 


(6) 


and hence 


h“ - 1 - 


yx'^dP/dy* )i 

y 0 ' 2 (dp/8y')o 


(7) 


It is to be noted that the structure of these equations, and all sub- 
sequent ones, is such that the quantities Pb> I > and which are pre- 
sumed to be given do not appear as independent variables. Instead, they 
are related parametrically to the shape and drag of the optimum profile 
through the parameter y’ . For example, equation ( 7 ) gives 
H (yx'> Yo’> M bo)> 8X15 equation 4 gives Pb(yi'> Yo'? Moo)> hence by select- 
ing arbitrary values of yx ’ the function H(Pb> Vo* > ^ 00 ) can be deter- 
mined. Moreover, as will be seen later, I depends on YoS and H, 
so thqt the parametric structure of the equations ultimately yields the 
desired function H(Pb, I, Mj, 


The differential equation (3a) and the appropriate boundary condi- 
tions y(o) = 0, y(s) =. t/2, and y(s + l) = \,j2. determine 
function of y by a single quadrature., 


x as a 
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/dp/dy 1 fly 


/ dP/dy 1 dy _ 
/(t/2)* 1 - y^ 


J 

*/X (s + l) + ~~ 


upstream 

surface 


downstream 


t /2 ** ^/ 2 ) n " ^ surface 


The alg ebraic sign to be selected for each radical is determined by the 
requirement that x increases as y increases on the upstream surface, 
and as y decreases on the downstream surface. The above relation is 
more conveniently expressed in terms of the dimensionless quantities 
Y = y/(t/2), X = x/s, and L = l/s. After eliminating by evaluating 

the above integral at x = s, there is obtained for the dimensionless 
shape 

f S(Y) on upstream surface 1 


' L + i(Y) on downstream surface 

where §(Y) is defined as the quotient of two integrals 1 


5 00 = 


i 


5p/5y t dY/V i - y* 


The definite integral which normalizes |(Y), such that £ (l) = 1, will 
for sake of brevity be designated by j n . Thus 


Jn = 


1 - HP- 


Actually, for completeness the function |(Y) should be written as 
£(Y; J^a, n, y 0 ’) because it depends on the three quantities M^, n, and 
y 0 * as well as the variable Y. For brevity, though, it is written 
simply as £(Y). 


Since x is a double-valued function of y over the chord length, * 
£(Y) is also a double-valued function of Y. A sketch of a typical 
curve of Y versus £ is shown in figure 2. For. a given dimensionless 
ordinate Y, one of the two values of § represents the dimensionless 
chordwise distance from the leading edge to a point on the upstream 
surface, whereas the other value represents (apart from an additive 
term L) the dimensionless chordwise distance from the leading edge to 
the point on the downstream surface which has the same ordinate Y. It 
is to be understood that in determining X from equation (8) the appro- 
priate value of £(y) must be used for each surface; thus, for a given Y 

^■Since the integrands are singular at the point Y = 1, numerical compu- 
tations of the function £(Y) must allow for this singularity. A 
simple method of doing this is outlined in appendix A. 
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the value of I (Y) appropriate to the downstream' surface always is 
greater than for the upstream surface. 


From the foregoing and the fact that £(y) is independent of Pb 
and a, it follows that a curve of Y versus £(y), such as is illus- 
trated in figure 2, determines the curved portions of an infinite number 
of optimum profiles, all having the same values of y 0 ’ , M^, and n. It 
maybe noted that a (equation (5)) essentially determines the length 
of straight midsection which is to be placed between the two curved 
portions after separating them at the point where Y = 1 (point AB in 
fig. 2)j whereas Pb (equations (4) and ( 7 )) essentially determines H, 
the value of Y beyond which the downstream portion of the curve is not 
used in a given- case. It is noted that although the chord is fixed, 
the value of £ corresponding to the trailing edge is not. This is 
because £ = x/s changes whenever s changes. 

As may be deduced from equation (7), yi T determines H for a 
given n, M,,,, and y 0 '. Moreover, y x ' determines Pb from equation (4). 
Hence, for any given value of base pressure the point on the downstream 
surface which corresponds to the trailing-edge position can be indicated 
on each |(Y) curve. (See fig. 2 where the point corresponding to zero 
base pressure is indicated.) 

Solution for L, ?, t/c, and s/c .- Turning now to the determination 
of the optimum length of straight midsection Z, one sees from equa- 
tion ( 5 ) that such a determination will also give T. Since Z is a 
function of yo* and the given quantities n, a, Pb, and M*,, this 
enables y 0 ', the quantity used as a parameter in the present analysis, 
to be related to the quantity I, which is a more convenient one to use 
if I is the actual quantity given. Starting with the definition of I 
the following equations result; 

I s l/(t/2) n “ ff 

= - / C ^Y n dX = -f /3 [1 - (1 - Y n )] dX 

c J o c J o 


s 
c 

or, by using the relation 

dX _ 1_ /BP/By 1 
dY = J n V 1 - Y n ; 

which follows from equations (8), ( 9 ), and (10), there results 


~ -Z? (1 - 1 “) £ ay 


dX 


dY 


1 - I 
s/c 



Bf 

5y» 


(1 - Y n ) 


dY 


( 11 ) 
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For convenience the right-hand member of this equation is defined 
as rj(H). 


7(H) 


s i_ r* /ST 

«3n Jq v dy' 


(1 - Y n ) dY 


( 12 ) 


From equations ( 8 ) and (5), 


7 = L + | (H) 
s 


Combining this with the above relation between I and s/c gives 


and 


s n - a 

c n£(H) - cft|(H) 


•(13) 


T = S(H) - 7(H) 

*(=)-£ 1(H) 


(U) 


Recalling that H is determined by Pb for a given !(y) curve (given 
M eo n > Yo’)j one sees that equation (llf) determines I as a function 
of n, o, Pb, and y 0 f . A convenient determination of I, of course, 

can only be made if the function t](H) in addition to 5 (y) has been 
computed. The function t](h), which for completeness should be written 
as q(H; M^, n, y Q '), is somewhat easier to compute than £(y) since it 
is not singular at Y = 1. Attention is called to the fact that all the 
above integrals with limits ranging from Y = 0 to Y = H, as in equa- 
tion ( 12 ), for example, really correspond to integration over both 
curved surfaces, first from Y = 0 to Y = 1, and then from Y=1 to Y = H. 


With the position of maximum thickness determined by equation (13), 
the maximum thickness ratio can be determined in terms of the surface 
slope at the leading edge. 


'dY\ _ 2s /dy\ 
d V o ~ t [dxjo 


or, from equation (ll). 


t 

c 



Yo 1 

(dY/dX ) 0 


in an alternate form 


t_ 

c 



} ( 15 ) 


( 15 a) 


Calculation of pressure drag of an optimum -profile .- In reference 4 
it was shown that for linearized supersonic airfoil theory the pressure 
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drag coefficient of an optimum profile was a simple algebraic function 
of certain quantities such as H and s/c. Since these quantities are 
known once the shape of the optimum profile is determined, a separate 
integration is not required in order to calculate the pressure drag. 
Fortunately, a similar algebraic relation can also be deyeloped for the 
present case. In so doing, integration by parts is employed starting 
with the defining equation for pressure drag. 


f ea - £ /2 pay - Pb § - Pi| yap - Pbf (16) 

In these equations, and subsequent ones, the integration is carried out 
only over the two curved portions since the straight midsection can con- 
tribute no drag. From equation (3a)> 


dP = 


dP 

d(l/y*) 


F) 


d Mr = -y- 


dp_ 

dy» 


-x [(t/2) n - /] a(§£) 


d 



hence, substituting into equation (l6) and again integrating by parts 
gives 


c 

2 



- ^ + v x 




[(t/2) n - (n + ljy^dx - 


x ^ b+i ^ t / 2 ^ n - ( n + 1 )y n 3 dx 


The first bracketed term on the right side vanishes by virtue of equa- 
tion (4). The remaining Integrals can be simplified by noting that 




r.C 


y 31 


dx, I =crj, and x(t/2) n = y 0 '‘ 


(g-i 


There results 


c d = 2y 0 ' 2 ^^) [(n + 1 - cr)T - l] (17) 

This equation enables the pressure drag to be readily calculated if the 
base pressure is given, since P^ determines H for a given y Q ’ and 
Moo, and H determines T in accordance with equation (lit-). Thus, 
equation (17) involves the base drag implicitly, but not explicitly. 

For the special case of linearized supersonic flow, 


dP 

dy 1 


= 2//M^ 


1 


r> 




and the above equation for c d can be shown to reduce to the correspond- 
ing equation for pressure drag developed in reference it-. 
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Closed-Form Solution for the Special Case of Given 
Cross-Section Area n = 1, a - 0 . 




When the cross-section area of a profile is prescribed (n = 1, 

0 = 0), corresponding to a given wing volume, torsional stiffness, or 
torsional strength of a thin-skin structure, then the differential 
equation (3) can be integrated immediately with respect to x to yield 
a solution in closed form for the airfoil shape. There results 


-Xx = P + y’ tr— - + constant 
ay* 

The constants can be eliminated by evaluating this expression at 
and x = s to obtain 


(18) 


x = 0 


f - x - «*> - 1 - p : p(0 > T 

8 P 0 + y 0 ’(dP/V) 0 - P(0) 


(19) 


Here P(o) is the pressure coefficient at y* = 0, and P Q is the 
pressure coefficient at x = 0. For practical purposes P(o) usually 
can be taken as zero, although strictly speaking it should be regarded 
as a small quantity compared to P Q + y 0 * (<dp/Sy' ) Q . The parametric 
equations for Y and H in terms of y* are the same as before, only 
with n ss 1. 


Y _ 1 _ y* g (8p/8y» ) 

y 0 f 2 (dP/V)o 


( 20 ) 


h = 1 - ( 21 ) 

y 0 ' 2 (8p/8y') 0 

The general equation for the base pressure coefficient does not involve n, 
and hence is the same as before. 

Pb = + V « 

The constant in equation (l8) can be evaluated at x = c i ns tead 
of at x * s. Combining such an evaluation with the above equation for P^ 
yields the alternate expression 


x u 1 . p + y(ap/ay) - Pb (22) 

c P 0 + y 0 f (ap/V ) 0 -Pb 

which involves P^ instead of P(0). The equations for t/c, s/c, and eg 
can be derived easily from the preceding equations. Omitting algebraic 
details, the following results are obtained: 
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By evaluating equation (l8) at x = c, and combining -with equations (3b) 
and (4), 


t - 2yo tg (5P/8y t ) 0 
“ " P 0 + y 0 *(dp/dy») 0 - Pb 

by evaluating equation (22) at x = s , 

B , P(0) - Pb 

c ' p 0 + y 0 , (^ p /^y , )o - ?b 

and by substituting n = 1 and a = 0 into equation (17), 

*• - 2 V 2 (|r) o <* - 1) 


(23) 


(24) 


(25) 


Solution for the Special Case of Given Thickness Ratio 

n = 00, a = finite 


In reference 4 it was shown that the limiting values n = 00 , 
a = finite, represent the auxiliary condition of a given airfoil thick- 
ness ratio. The mathematical simplification inherent in the use of 
approximate theories such as linearized flow enables the solution for 
a given thickness ratio to be obtained directly by passing the general 
solution to the limit as n ->«. For shock-expansion theory, though, 
a general’ solution in closed explicit form cannot be obtained, and 
recourse to the alternate method indicated in reference 4 is required. 
This alternate method deals directly with the appropriate differential 
equation, which, for the case of given airfoil thickness, becomes simply 

yi 2 = constant (26) 

oy* 

which is satisfied by a profile composed of : any number of straight seg- 
ments. As shown in reference 4, the constant in the above equation does 
not change over the entire chord, with the result that the upper half 
of the profile forward of the 'trailing edge is composed of two straight 
lines, one extending from x=0 =ytox=s, y = t/2, and the other 
extending from x = s, y = t/2 .to the trailing edge x = c, y = h/2. 

The slope is discontinuous at the point where y = t/2.- To obtain a 
solution using any given airfoil theory, it is necessary to satisfy the 
differential equation (26), the end condition (4), and the boundary 
condition of a fixed thickness ratio. 

* 

Equation (17) for c<j becomes indeterminate as n—> 00 because 
T*H>0. For this case, however, the shape is known and the pressure drag 
can be determined from simple physical considerations: 

(=a) glTC n t/c - * t p o - rR 1 - H > - p b Hl 
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APPLICATION OF ANALYSIS 
Flow at Infinite Mach Number 


Prior to considering shock-expansion theory, the relatively simple 
case of flow at infinite Mach number over slender airfoils with small 
surface curvature will be considered. For such conditions the pressure 
coefficient on a surface facing upstream is proportional to the square 
of the local surface slope, that is, P = Gy' 2 . Since the pressure coef- 
ficient is zero on any surface facing downstream it follows from physical 
considerations that H = 1. Equation (4) is satisfied by requiring that 
yj_’ = 0. By consideration of the differential equation (3a) as special- 
ized to the present case it follows that 



= (2c) 


1/3 r 


L 


1 t 


n" 


1/6 


(1 - Y n ) 


l/e 


By substitution into equation (9), and employment of gamma functions to 
evaluate the integral in the denominator, there results 



( 27 ) 


The function Y versus |(Y) is plotted in figure 3 for n = 1, 2, 3, and «. 
The infinite value of n corresponds to the auxiliary condition of a 
given thickness ratio, and the optimum profile in this case is a wedge, 

since and. S(Y)->Y as n->». It is seen that there is 

little difference between the three curves for finite n. 


The other characteristic function needed for the complete determina- 
tion of an optimum profile Is t)(h). By substituting H = 1 in equa- 
tions (12), ( 13 ), (14), and (17), and employing equations (5), (10), and 
(28), the following expressions are obtained: 


3 

(28) 

s _ (n - cr)(2n + 3) 

(29) 

c n(2n + 3 - 2a) 


I 

(30) 

2n + 3 - 2a 

1 3ct 

c n(2n + 3 - 2a) 

( 31 ) 



Ik 
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t_ 

c 



2 

c 


(3 + 2n - 2a) 
3 


I. 


1 

n-q 



(32) 

(33) 

m 


If desired, this last equation for the pressure drag coefficient can be 
■written in terns of I instead of t/c, inasmuch as I is related to 
the thickness ratio through equation (32). It should be noted that the 
above equations relate in closed form all pertinent properties of the 
optimum profile to the given quantities I, n, and 0 . Examples of opti- 
mum profiles determined with the aid of these equations are presented 
subsequently. 


Shock-Expansion Theory 


When the oblique shock-wave and Prandtl-Meyer equations are combined 
to calculate the pressure on an airfoil surface in supersonic flow, the 
resulting equations for P are quite involved. The appropriate equa- 
tion for 8p/8y* , however, may be obtained by starting with the local 
differential relation 


1 8p = 2 . 

|pV 2 86 Jk* - 1 


(35) 


This point relation is formally the same as the corresponding relation 
applied throughout an entire flow field in linearized supersonic airfoil 
theory. The partial derivative 8p/8b is * taken with and 5 0 held 
constant. Expressing equation (35) in terms of the pressure coefficient 
and free-stream conditions, 

&P = pV g 2 . 

& ~ PoqVoo 2 VM 2 -1 


or, since 


8p 

8y' 


5p 

85 


cos 2 5 


8P ■ pV 2 2 cos 2 S 

dy' 8^0?" -s/m 2 - 1 


there results 
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This equation appears fairly simple, but is not in a form which, can be 
computed readily with the aid of existing tables where quantities such 
as local pressure ratio p/pt and total-pressure ratio across a shock 
wave Pt/Pt are tabulated. Thus, a more convenient form for calcula- 
ting purposes is 


. g— [Ml / M » - 1 /pfelfELWe 
V m w 2 - 1 LmJ v m 2 - 1 (p M /ptJ \Pt J 


(36) 


From this equation a numerical value of Bp/cty* can be determined from 
tabulated oblique-shock and expansion characteristics once y 0 ’, M^,, 
and the local slope y’ are specified. 


The functions £(Y) and tj(H) have been calculated for shock-, 
expansion theory by substituting equation ( 36 ) into equations ( 9 ) and (12), 
respectively, and then performing the indicated integration graphically 
by the method outlined in appendix A. In this process other useful 
quantities axe calculated such as (Sp/5y* ) Q and j n . The results are 
presented in figures k , 6, and 7. In figure 4 the quantity 

(dP/dy* ) 0 /(2 //m„ 2 - 1) , which is equal to (dp/dy f ) 0 /[(d?/dy* ) 0 ] g 

is plotted as a function of 8 0 for various values of M B . Similarly, 
<3n/(Jn)s =o° ls P lotted in figure 5 . It is to be noted 


0 ° 


where 




r 1 dY 

- 2. Jo J l - 


2kn 


</ M»2 - 1 


k 


n 


2 for n = 1 
rt/2 for n = 2 
1.4023. ..for n = 3 


(37) 


Curves of H 31 versus 5o for various values of base pressure are pre- 
sented in figure 6, from which it is apparent that the dimensionless 
trailing-edge thickness increases if either the airfoil thickness 
increases (5 0 increases), or if the base pressure increases. In figure 7 
the functions |(Y) and tj(H) for various and n are presented 

plotted in the form Y versus §, and H versus q. The curves of Y 
versus | determine the shape of the optimum profile, while the curves 
of H versus q are useful in determining I, s/c, and t/c. Many of the 
curves of Y versus I have been terminated at the point (indicated by 
small circle) corresponding to zero base pressure. 


EXAMPLES AND DISCUSSION 


In order to determine an optimum profile it is necessary, of course, 
to know the base pressure. Experiments have shown that base pressure in 
two-dimensional flow depends principally on the Mach number, type of 
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boundary- layer flow, and the boundary- layer thickness at the base. 

(See reference 6.) Average experimental values are shown in figure 8 
for both laminar and turbulent flow plotted as a function of the parame- 
ters proportional to the ratio of boundary- layer thickness to trailing- 
edge thickness. Step-by-step details of the method of determining an 
opt imum profile by combining experimental base pressure data with the 
curves of figures 1). to 7 are given in appendix B. 

In figure 9 examples of optimum profiles determined by the theory 
of the present report are shown together with corresponding profiles 
determined by linearized theory (reference J+-). For each of the various 
auxiliary conditions the particular value of I selected for these 
examples is equal to that for a circular-arc biconvex airfoil of thick- 
ness ratio t C a/ c = 0.0 6. Since ..... — 

2 2n (n!) 2 (t C a/2) n - a 
Ica " (2 n+l)l 

(see reference k), it follows that with c = 1 the optimum profiles in 
figure 9 correspond to the value 

„ 2 2n (^) 2 (0.03) n “ g 
1 ” (2n + .1)1 

The auxiliary condition for n * 00 (fig. 9(d)) corresponds to a given 
maximum airfoil thickness of 0.06c. As indicated in figure 9 (and also 
in subsequent figures of this report), the base pressure for M#* = 1.5 
and M w = 3.0 corresponds to turbulent boundary-layer flow at a Reynolds 
number of 10 Since h is involved in the abscissa of figure 8, due 
allowance is made for the variation of base pressure with trailing-edge 
thickness. Because base pressure data are not available as yet for 
M„ = 8, a constant value has been assumed (pb = 0.1 p*) which is believed 
to be reasonable for a moderately thick trailing edge (h/c r* 0.05 or more) 
but probably greatly overestimates the base drag for a thin trailing edge 
(h/c ~ 0.01 or less). For 1^=00, it is not necessary to know the base 
pressure since the optimum profile at this limit is independent of Pb* 

In figure 10 examples are. shown for various values of I with 
Mca = 3, n = 3, and a = 0. Instead of specifying the value of I in each 
case, the thickness ratio of a structurally equivalent circular-arc 
biconvex airfoil is specified, as the significance of this latter value 
is easier to visualize. As would.be expected, there is no appreciable 
difference between the profiles determined by linearized and shock- 
expansion theory when the value of I is small (tea / c = 0.02, for. 
example), although differences are evident for larger values of I 
(t ca / c = 0.04 and larger). 

For each of the examples shown in figures 9 and 10, it is to be 
■ noted that in comparison to the profile determined by linearized theory, , 
the corresponding optimum profile determined by shock-expansion theory 
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has a smaller slope over the portion of surface facing upstream, a 
position of maximum thickness farther aft, and a greater slope over the 
portion of surface facing downstream. This difference which increases 
with increasing Mach number is to be expected, as indicated in refer- 
ence 4, because the linearized theory overestimates the suction forces 
and underestimates the positive pressure forces. 

As regards drag coefficient, it is evident that drag calculations 
based on linearized theory cannot be used at hypersonic Mach numbers 
since the computed coefficient approaches zero as the Mach number 
increases. Also, it is to be remembered that linearized theory is con- 
siderably less accurate in predicting the drag of blunt-trailing-edge 
airfoils than of sharp-tralling-edge airfoils, since the Busemann second- 
order terms for the upstream and downstream surfaces do not cancel as 
they do when the trailing edge is sharp. 

If shock- expans ion theory is used to calculate the drag of the pro- 
file determined by linearized theory, the resulting value is only slightly 
greater than the drag of the same profile determined by shocpc-expansion 
theory. In order to put this idea on a more firm quantitative basis, 
the case of inf inite Mach number can be considered, as the differences 
between profiles determined by' linear and nonlinear theory are the 
greatest at this limit. (See examples in fig. 9.) By use of the expres- 
sion P = Cy’ 2 to calculate the drag of the profile determined by linear 
theory, and dividing by the drag calculated correspondingly for the pro- 
file determined by nonlinear theory, the following expression results: 

'1.08 for n=l, 0=0 

1.06 for n=2,a=0 

1.07 for n=2,cr=l . 

1.04 for n=3,cr=0 

1.05 for n=3,a=l. 

^1.00 for n=co,cx finite 

(38) 





n+2-cr 


l6(n+l) 

_(n+2)(2n+3-2a)_ 



2(2n+3-2g) 

3(n+2-or) 


a 

n-a 


2n(2n+3)r 


!*t) 


G) r C n ) 


It is seen that at infinite Mach number the actual pressure drag of an 
optimum profile whose shape is determined by linearized theory does not 
exceed the drag of the true optimum profile by more than about 8 percent. 

If consideration Is given to the consistent differences noted earlier 
between the shapes of optimum profiles determined by linear and nonlinear 
theory, it is evident that linearized theory can be used with good 
accuracy to determine the optimum profile at any supersonic Mach number 
up to infinity. 2 Even without considering the consistent difference noted 
above, the profile determined by linearized theory is sufficiently accurate 
for most engineering purposes. Under less general condit ions a similar 

2 As indicated in reference 4, however, the linearized theory does not 
yield a reasonably accurate profile at the low supersonic Mach numbers 
near or below shock detachment. 
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result also has been found in the recent investigation of Klunker and 
Harder (reference 7) which appeared while the present report was being 
prepared. The shape of some of the optimum profiles determined in 
reference J, however, does not agree with the shape of analogous profiles 
in this report. For' example, it is indicated in reference 7 that the 
profile of least drag for t/c = 0.06 has a : sharp trailing edge at all 
Mach numbers below about 6, whereas the corresponding profiles shown in 
figure 9(3) indicate appreciable trailing-edge thickness even at Mach 
numbers of 1.5 and 3. This discrepancy is attributed to the arbitrary 
base pressure curve assumed in reference 7 which does not con-espond to 
measured data for thin trailing edges. — - 

From an engineering viewpoint it is desirable to know how much lower 
the drag of an optimum profile is than that of a sharp-trai ling-edge 
profile, and also how much the optimum profile can be altered without 
significantly increasing the drag. In order to provide a basis of com- 
parison, the zero-lift pressure drag of a family of sharp-trailing-edge 
circular-arc biconvex airfoils of various thickness ratios has been 
calculated by shock-expansion theory for the Mach number range between 
1.5 and 8. The results are shown in figure 11. Thus, for any profile 
the drag of a structurally equivalent (same value of I) circular-arc 
biconvex profile can be determined readily from the curves In figure 11 
by simply calculating t C a from the equation 

I * I ca “ ^(n !) 2 (t ca /2)“-<7(2n + l)l 

Computations of drag have been made for a family of "Bemioptimum" 
profiles having arbitrarily selected values of trailing-edge bluntness H, 
a shape forward of the trailing edge that yields minimum foredrag for each 
particular H, and the same value of I as a circular-arc biconvex pro- 
file of thickness ratio tea. These calculations have been carried out 
for t ca /c = 0.02, 0.0k, 0.06, 0.08, and 0.10 at Mach numbers of 1.5# 3> 
and 8, and for various combinations of n and 0. As In previous 
examples, the base drag in each case was determined from the curves of 
figure 8 for turbulent-boundary- layer flow at Re = 10 7 . The results are 
shown in figure 12 plotted in the form of a drag ratio versus H. Each 
curve corresponds - to a constant value of I, and is identified by the 
thickness ratio (t ca /c) of a circular-arc biconvex profile having the 
same value for I. In order to maintain a constant value of I, the 
actual thickness ratios (t/fc) of the semioptimum profiles change somewhat 
as H varies between 0 and 1 (the ratio t/t ca lies between about 0.90 
and 1.05 for the case n = 1, 0 = 0, between about 0.71 and 0.8k for 
n = 2, c=l, and between about 0.97 and I.08 for n = 3> 0 = 0) . For 
each curve in figure 12 the semioptimum profile having the minimum drag 
coincides with the optimum profile determined from the curves of fig- 
ures 5 to 8. The ordinate of each minimum point indicates the relative 
drag of the optimum compared to a structurally equivalent circular-arc 
biconvex profile, while the rise on each side of the minimum indicates 
the drag penalty resulting from the use of too much or too little 
trailing-edge thickness. It may be noted that some of the curves do not 
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cover the complete range of values of H. In all such cases, however, 
sufficient calculations were made so that the min imum point was bracketed. 

Perhaps the most significant feature apparent from the curves of 
figure 12 is the large increase in drag that results for thin airfoils 
when a trailing-edge thickness considerably greater than the opt imum is 
employed at moderate supersonic Mach numbers . The pressure drag of full- 
blunt profiles at Mach numbers of 1.5 and 3 for the case t C a/ c = 0.02 is 
several times the pressure drag of the optimum, the exact factor varying 
between about 2.8 and 3.9 depending on M^ n, and a. On the other 
hand, for thicker profiles (t C a/ c ^ 0 . 08 ) the pressure drag is much less 
sensitive to variations in trailing-edge thickness about the optimum, 
and the use of a full-blunt profile instead of the optimum would result 
in a much smaller percent drag penalty. Also, for a given value of t ca /c 
it is evident from figure 12 that the pressure drag becomes less sensitive 
to variations in trailing-edge thickness from the optimum as the Mach num- 
ber is increased into the hypersonic regime. 

As would be expected, the semioptimum sharp-trailing-edge profiles 
(H = 0 in fig. 12) have somewhat less pressure drag than a structurally 
equivalent circular-arc biconvex profile. The observed difference in 
drag for thin airfoils at moderate Mach numbers is negligible for the 
case of n = 1, a - 0, since the optimum sharp-trailing-edge profile for 
these conditions is very close to a circular-arc biconvex profile (if 
linearized theory were employed the optimum sharp-trailing-edge profile 
for n = 1, a = 0 would be a circular -arc biconvex profile). The 
corresponding difference in drag for the case n = 2, a = 1, however, is 
significant since the optimum profile in this case has a midsection of 
constant thickness, and hence is of considerably different shape, as 
well as being considerably thinner than a structurally equivalent circular- 
arc biconvex profile. 

Each of the curves for M. = 1.5 and 3 (figs. 12(a) and 12(b)) show 
a minimum at same finite value of H, but the curves for t ca /c = 0.02 at 
M = 8 (fig. 12(c)) do not. If all other parameters were constant, this 
trend would not be expected inasmuch as the optimum trailing-edge blunt- 
ness for a given thickness ratio generally increases as the Mach number 
is increased into the hypersonic range. (See reference 4.) The unexpected 
trend is observed in the present examples because the variation of base 
pressure with trailing-edge thickness is considered at M = 1.5 and 3 
where experimental measurements are available, but? it is not considered 
at M» - 8 where, in the absence of experimental data, a constant base 
pressure was arbitrarily assumed (one-tenth of the free-stream pressure, 
irrespective of trailing-edge thickness). It is expected that if base 
pressure measurements were made at M^ = 8, they would show a dependence 
on trailing-edge thickness just as at the lower Mach numbers. Conse- 
quently, it is believed that the actual curves for the thinner airfoils 
at U m = 8 will be greatly different than shown in figure 12(c), although 
the curves for the thicker airfoils are not expected to be significantly 
different. If a constant base pressure corresponding to measured values 
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on thick trailing edges were assumed at = 1.5 and. 3> the curves 
for t C a/c = 0.02 would rise starting from H = 0 even'more steeply 
than at = 8. This illustrates the necessity of considering the 
dependence of base pressure on trailing-edge thickness in such an 
analysis. 3 

The drag of the various optimum profiles (minimum points in fig. 12) 
is seen to be less than the drag of a structurally equivalent circular- 
arc biconvex profile by amounts varying between about 1 percent and 
53 percent, depending on the values of n, cr, M*, and t ca /c. The largest 
drag reduction occurs for the case of n = 2, a = 1, just as indicated 
by the linearized theory of reference it-. Likewise, the drag reduction 
generally increases as the Mach number or the thickness ratio is 
increased, as predicted by linearized theory. 


CONCLUSIONS 


1. For a given Mach number and structural requirement, the shape 
and drag of the profile having the least possible pressure drag at zero 
lift, as computed from shock-expansion theory, can be determined readily 
provided the base pressure is known from experiments and provided curves 
of certain characteristic functions are available. (These functions 
are |(Y, 6 0 , n), t)(H, 6 0 , n), H(& 0 > Pb/P«>, n), j n (&o)> and dp/dy' (B 0 ), 
examples of which are presented in figures 4 to 7.) 

2. A comparison of profiles determined by shock-expansion theory 
and linearized theory indicates that the linearized theory may be Used 
with reasonable accuracy at Mach numbers up to infinity to determine the 
shape of the optimum profile, although it can be used only at moderate 
supersonic Mach numbers to determine the drag. 

3. Considerable deviations in profile shape from the theoretical 
optimum can be made without increasing the drag excessively provided the 
Mach number is high, or the airfoil thickness ratio is relatively large. 
Large drag penalties result, hpwever, if a trailing-edge thickness 
appreciably greater than the optimum Is employed on a thin airfoil at 
moderate supersonic Mach numbers. 

4. It is necessary to consider the experimentally observed depend- 
ence of base pressure on trailing-edge thickness when calculating the 
optimum-profile shape and drag of a thin airfoil. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Moffett Field, Calif., June 25 j 1952 

3 If the boundary layer were laminar, the effect on the optimum shape of 
the dependence of base pressure on trailing-edge thickness would be 
even greater. 
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APPENDIX A 

METHOD OF CALCULATING |(y) 


The singularity at Y = 1 of the integrands in equation ( 9 ) causes 
difficulty when directly evaluating such integrals numerically or graphi- 
cally. This difficulty, however, can he circumvented by transforming 
from Y as the integration variable to a new function fn(Y) defined as 
follows : 


fn(Y) 



dY 

Vl - Y D 


(Al) 


With this transformation the equation for £(y) becomes 


In this 
finite, 
ence L, 



(A2) 


equation the integrands and the ranges of integration are all 
The co ns tant k n is the same quantity as that used in refer - 
namely. 


f 2 


kn = fn(l) 


) «/2 
1.4023 
. 1 


for n = 1 
for n = 2 
for n = 3 

for n = 00 




(A3) 


The functions fn(Y), apart from an additive constant, are likewise the 
same characteristic functions as appeared throughout the analysis of 
reference 4 when linearized supersonic airfoil theory was employed for P. 
Thus, 

f 2(1 - Vl - Y n ) for n = 1 1 


fnW 


= 


sin" 1 Y 

1 

1.U023 - 3" 4 F(k,9) 


for n = 2 
for n = 3 


(AL) 


Y 


for n = co 


j 
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where F(k, 9) is the incomplete elliptic integral of the first kind of 
modulus k = sin 75° = 0 . 9659 * 8^3 amplitude .. 


9 = 


cos -1 


1 + y) 

(v r 3 + 1 - Y) 


With the transformation to ^n> the integrals in equation (A2) are 
evaluated by first selecting a number of values of y 1 ranging from y 0 ’ 
to large negative yalues. For each y' the_ ordinate Y is computed 
from equation (6), f n from equation (A4), and 6p/8y* from the particu- 
lar airfoil theory. A plot is then made of cSP/By’ versus fn in order 
to evaluate the integrals determining i(Y). 
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APPENDIX B 

DETAII5 OF METHOD OF DETERMINING THE SHAPE AND DRAG OF 
AN OPTIMUM PROFILE BY SHOCK-EXPANSION METHOD 


In the shock-expansion equations the leading-edge deflection 
angle 8 0 is a more convenient parameter to use than the given value 
of I, hence the steps outlined below involve an iterative procedure. 

(1) Assume values of 5 0 and p^/p^ 

(2) Read value of H n from figure 6 and compute H 

(3) Read values of £(H) and T)(h) from figure 7; compute 1 from 

equation ( Ih-) and s/c from equation ( 13 ) 

(4) Read value of (dp/dy’ )o from figure 4, j n /0n)& o = o° 

from figure 5; compute (Jn)g 3 qo from equation (37 ) > 
t/c from equation ( 15 a) 

(5) Compute I = I (t/2) n-G 

By comparison of the computed value of I with the given value, a new 
value of 5 0 can be estimated. Also, from the computed value of h/c, 
the experimental base pressure curves in figure 8 yield a new value of 
Pb/P w - By repetition of the above steps until the final computed value 
of I Is equal to the given value, and the final computed value of h/c 
corresponds to the final base pressure assumed, all characteristics 
(8 0 , t/c, s/c, H, |(Y), (dp/dy* ) Q , and 7) of the optimum profile are 
determined. The pressure drag is then calculated from equation (17). 
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y^tanS 0 

yi=tan$, A B 



Figure /.- Sketch illustrating various symbols used in the analysis. 



Dimensionless chordwise distance , £(Y) 


Figure 2. - Sketch of a typical curve for determining optimum-airfoil 

profiles. 
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Figure 3.- Characteristic curves applicable to the simplified equation for flow at infinite 

Mach number. 
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Figure 4.- Curves based on shock-expansion theory showing (f^/(2/jMSr ) 
as a function of leading-edge deflection angle for various Mach numbers. 
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Figure 5. — Curves based on shock-expansion theory showing 
as a function of leading-edge deflection angle for various 
numbers 
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Figure 6. — Curves based on shock-expansion theory showing dimension- 
less traiiing-edge thickness as a function of base pressure and 
leading-edge deflection angle. 
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Figure 6 . — Concluded. 
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Figure 7 -Continued. 
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Figure Z —Continued. 
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c/[h(Re) ,/z ] 



c/[h(Re) ,/s ] 

(b) Turbulent flow. 

Figure 8 .- Average experimental values of base pressure from reference 6 . 



shock expansion theory 
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Figure 9 Examples of optimum profiles f auxiliary integral equal to that of 
a circular-arc biconvex airfoil of thickness ratio 0.06. 
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(b) n* 2,o-* I 
Figure 9 Continued. 
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(c) n* 3 , «r a 0 

Figure 9 - Continued 
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Shock -expansion theory 
Linearized theory 
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Figure 9- Concluded. 
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Shock- expansion theory 
linearized theory 






Figure 10 -Examples of optimum profiles for values of auxiliary integral equal 
to that of a circular- arc biconvex airfoil of maximum thickness tea-, n=3, 
cr-O-, M ao = 3 i base pressure corresponding to turbulent flow at Re* 1 0 7 . 
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Figure 1 1 .-Pressure drag of circular-arc biconvex airfoils determined by shock- 

expansion theory. 



(a)M p *i.5, p b corresponding to turbulent flow at Re*IO ? 

Figure 12.- Effect of trailing -edge bluntness on the relative pressure drag 
of semioptimum profiles compared to a structurally equivalent circular- 
arc biconvex profile. 
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(b) Mgf 3, p b corresponding to turbulent flow at Re * 1 0 7 . 



(clMa'S, PfO.lp^t assumed). 


Figure 12. — Concluded. 
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